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ABSTRACT 

We introduce a set of two-parameter models for the dark energy equation 
of state (EOS) w{z) to investigate time-varying dark energy. The models are 
classified into two types according to their boundary behaviors at the redshift 
z = (0, oo) and their local extremum properties. A joint analysis based on four 
observations (SNe + BAO + CMB + Ho) is carried out to constrain all the 
models. It is shown that all models get almost the same Xmin — 469 and the 
cosmological parameters {Qm, h,flbh'^) with the best-fit results (0.28,0.70,2.24), 
although the constraint results on two parameters {wo,Wi) and the allowed re- 
gions for the EOS w{z) are sensitive to different models and a given extra model 
parameter. For three of Type I models which have similar functional behaviors 
with the so-called CPL model, the constrained two parameters Wq and Wi have 
negative correlation and are compatible with the ones in CPL model, and the al- 
lowed regions of w{z) get a narrow node at z ~ 0.2. The best-fit results from the 
most stringent constraints in Model la give {wo,Wi) = (— 0.961q;2i, — 0.12^Q;gg) 
which may compare with the best-fit results {wo,Wi) = (—0.97^0^1, —O.lb^i'^l) 
in the CPL model. For four of Type II models which have logarithmic func- 
tion forms and an extremum point, the allowed regions of w{z) are found to be 
sensitive to different models and a given extra parameter. It is interesting to 
obtain two models in which two parameters Wq and Wi are strongly correlative 
and appropriately reduced to one parameter by a linear relation wi oc (1 -|- wq). 
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1. INTRODUCTION 

The distance determinations at cosmological scales, from standard candles such as Type 
la supernovae(SNe) (Perlmutter et al. 1999, Riess et al. 1998) as well as standard rulers 
such as the cosmic microwave background(CMB)(Hinshaw et al. 2009; Komatsu et al. 2009) 
and the baryon acoustic oscillation (B AO) (Eisenstein et al. 2005; Cole et al. 2005; Huetsi 
2006; Percival et al. 2007), have firmly established that the expansion of our universe in 
the time close to the present was accelerating. It is now believed that a mysterious negative 
pressure component in the universe, named dark energy, is responsible for this accelerated 
expansion, which contributes about 70% critical density of universe today. The underlying 
physics of dark energy remains obscure and most of relevant studies parameterize dark energy 
by its equation of state (EOS) parameter w = p/p. In observational data analysis, it is 
widely assumed that the dark energy EOS is independent on specific theoretical models and 
follows a certain predetermined evolutionary history. Common parameterizations include 
the constant the linear variation of the redshift w{z) — Wq -\- w^z, a convergent form 
w{a) — Wo + Wa{l — a) with a — 1/(1 + z) the scale factor, and the piecewise functions binned 
on the redshift. Despite the fact that the accumulated observational data have significantly 
improved the constraint on the dark energy EOS, our knowledge about the time evolution 
EOS w{z) is rather poor in considering its infinite dimensional parameter space. Thus, it is 
useful to study more general classes of models with the time-varying dark energy. Here we 
shall utilize the 397 SNe la provided by Hicken et al. (2009), 7- year WMAP observations, 
the BAO data given by SDSS DR7, together with the Hubble constant Hq (Riess et al. 2009) 
constrained from the SHOES team to investigate varieties of the parameterizations of dark 
energy EOS in the fiat cosmology. 

By far, all observed data are consistent with the ACDM cosmology, with dark energy in 
the form of Einstein's cosmological constant A or the vacuum energy, which is time indepen- 
dent and has an EOS w = —1. However, this model suffers from the difficulties related to 
the fine tuning and the coincidence problem (see e.g. Padmanabhan 2003). In last decades, 
several kinds of dynamical dark energy models have been proposed, such as holographic, 
quintessence, phantoms, K-essence, branes, Chaplygin gas, etc. In these situations, dark 
energy EOS changes with time and is generally parametrized as w{z) = wq -\- Wi f{z), with 
Wq the present value of dark energy EOS for f(0) = 0. Conveniently, it includes the case of 
the ACDM model {wq = —l,Wi = 0), and the constant EOS {wq = w,Wi = 0). 

The constant EOS is the simplest choice of dark energy parameterizations and a well 
approximation of the slowly rolling quintessence field (Huey et al. 1999; Wang et al. 2000). 
When determining whether or not the dark energy is a cosmological constant, which is the 
primary goal of the study of dark energy, the constant EOS is the most direct form. Thus 
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almost all cosmological observations data provide their constraint on the constant w. By 
far its value is found to be very close to w = —1 and the latest limit given by WMAP 
together with BAO and Hq is w = —1.10 ± 0.14 (Komatsu et al. 2010). For the time- 
varying dark energy EOS, the simplest parameterizations is the first-order Taylor expansion 
w^°'^{z) = wq + w' z (Huterer and Turner, 2001; Weller and Albrechet, 2002), where wq and 
w' are the value and the first-order derivative of dark energy EOS at present, respectively. 
Though this model was used in SNe la data analysis, its divergence behavior at the large 
redshift 2; — > oo is, in principle, disfavored and at the same time, prevents it from CMB 
data fitting. The so-called CPL parameterization w'^^'-"{z) — wq-\- w z z / {1 -\- z) , proposed by 
Chevallier and Polarski (2001) and Linder (2003), overcomes this defect and has capacities 
to describe a variety of scalar field theories. For recent low redshift z « 1, w^^^ is 
equivalent to the 1-order Taylor expansion w^""^ mentioned above, and for high redshift, its 
value approaches Wq + Wz- By now, the CPL parameterization has been adopted widely for 
all kinds of cosmological observations. 

There are other remarkable works on the parameterization of dark energy EOS. The 
logarithmic expansion in z was proposed to approximate scalar field EOS of the quintessence 
models with smoothly varying potential (Efstathiou 2000; Gerke and Efstathiou 2002). The 
CPL-extended parametrization was proposed through changing the power exponent of the 
denominator in the CPL model from 1 to 2 (Jassal, Bagla and Padmanabhan 2005). Under 
this model, dark energy evolves rapidly at high redshift and has the same EOS at the present 
epoch and the low redshift. Several attempts have also been made in the past to parameter- 
ize dark energy EOS w{z) with more than two parameters. The kink approach utilizes four 
observables to parameterize dark energy EOS: the present value w{z = 0), the early value 
w{z — )■ 00), the scale factor at the transition point at and the width of the transition A (Bas- 
sett et al. 2003, Corasaniti and Copeland 2003, Corasaniti et al. 2004). This model is very 
useful when studying the rapid transition from dccclacratcd to accelerated expansion. Using 
the gold sample of SNe la (Riess et al. 2004) and a matter density prior. It was compared 
for the fitting results from the standard Taylor expansions of w{z) and the Kink parame- 
terization (Bassett, Corasaniti and Kunz 2004). For the form of Taylor expansion, it was 
considered with four different expansion functions and three different expansion orders, high- 
est to second order. With the marginalised on-dimensional likelihoods for the parameters, it 
has been demonstrated that expansion order is more important than paramctrisation. After 
illustrating the 'maximised' limits on the redshift of dark energy EOS w{z), it was claimed 
that the standard two-parameter expansions artificially rule out the models associated with 
rapid evolution of dark energy. When focusing on the observational information at high red- 
shift, a variation of the Kink parameterization was adopted by introducing the mean value 
< w{z) > of dark energy EOS as an explicit parameter (Pogosian et al. 2005). Different from 
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the two-parameter models by fixing the value of the EOS today and at high redshift with a 
definite evolution process, three-parameter or four-parameter models are more flexible and 
allow for both rapid and slow variations of w{z). 

In our past work(Zhang et al. 2009; Zhang and Wu 2010), we have present the constraint 
on the time-varying dark energy from the splitting angle statistic of strong gravitational 
lenses, combined with SNe la and BAO observations. In this paper, by taking the SNe 
la data(Hicken et al. 2009), the baryonic acoustic oscillations(Eiscnstcin et al. 2007), the 
7- year WMAP observations(Komatsu et al. 2009), together with the Hubble constant Hq 
(Riess et al. 2009), we shall investigate a set of two-parameter models of dark energy EOS 
w{z). We will mainly highlight two issues which have not previously been illuminated. First, 
we pay attention to the detailed study on some interesting two-parameter models with two 
kinds of boundary behaviors. Second, we shall carefully demonstrate that there exist two 
kinds of strongly correlative dark energy parametrizations so that the two parameters are 
reduced to one parameter. Our paper is organized as follows: The cosmological observation 
data used in our analysis are presented in Sect. 2. In Sect. 3, we make a joint analysis on 
a set of two-parameter models proposed in our present paper. The conclusions are given in 
the last section. 



2. Relevant Cosmological Observation Data 

In this section we describe four sets of different cosmological observation data used in 
our parameterization analysis, namely SNe la, BAO, CMB and Hubble constant Hq, which 
are the currently preferred probes of the expansion history of the universe. Besides the well- 
controlled system uncertainties and high measurement precisions, their analysis formulas are 
quite simple and time-saving. This is important for our numerical calculations when con- 
sidering the relevant cosmological parameters simultaneously. The H{z) observations from 
Stern et al. (2010) have also been investigated and their additional constraining power is 
not notable, and we shall not include the H{z) sample here. As the dark energy affects the 
space-time geometry and the distance determination, the SNe la sample and BAO measure- 
ments provide us the distance ladder at low redshift z < 1.5 and are expected to distinguish 
the varieties of dark energy EOS w{z) in the later epoch, while the much higher character- 
istic redshift z ~ 1090 of CMB shall have somewhat "integral" effect and at the same time, 
give the limit on w{z) in the early universe. Combining above three observations makes 
it reasonable to parameterize the dark energy EOS w{z) by its two special limiting values 
w{z = 0) and w{z — >■ oo). On the other hand, including the measurement of Hubble constant 
Hq can help to break the parameter degeneracies in BAO and CMB observations. 
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2.1. Type la supernovae 

It is widely believed that SNe la has homogeneous intrinsic luminosity of peak magnitude 
and thus the SNe la data provide the precise distance measurements at the cosmic scale. 
The analysis of the distance modulus versus redshift relation of SNe la shows us the direct 
evidence of the existence of dark energy. In the last decades, many SNe la observations have 
been done and the total number of SNe la sample increases quickly(e.g. Riess et al. 2004; 
Wood-Vasey et al. 2007; Kessler et al. 2009). The SN la compilations are often consist 
of high-redshift {z ~ 0.5) data set and low-redshift (z ~ 0.05) sample at the same time. 
Remarkably, Kowalski et al. (2008) provided the Union data set with a compilation of 307 
SNe la discovered in different surveys. The heterogeneous nature of the data set have been 
reflected and all SNe la sample are analyzed with the same analysis procedure. Then by 
adding the CAF3 sample, Hicken et al. (2009) shown the Constitution set with 397 SNe la, 
of which all SNe la light curves are fitted by using the spectral-template-based fit method 
of Guy et al. (2005) (also known as SALT). In this paper, we use this constitution set to 
constrain the dark energy EOS. 

In the fiat universe, the Friedmann equation are given by 



H(z)/H„ = V^mCI + zf + (1 - J!„)/(z) 

/(.) = a.p(3y['il^d.') (1) 

with Hubble constant Hq — 100 h km s~^ Mpc~^. The influence of dark energy EOS w{z) 
is focused on the dark energy density 0,de{z) = (1 — ^M)f{z)- For theoretical calculations, 
the luminosity distance of SNe la is defined as 

In the SNe la observations, it is traditional to provide the distance moduli /j,, which is related 
to the above luminosity distance c/l by 

// = 51ogdi/Mpc + 25. (3) 

Due to the entire degeneracy of the Hubble constant Hq and the absolute magnitude of SNe 
la, the parameter Hq is irrelevant for the SNe la data analysis. Thus we use the following 
statistic (Perivolaropoulos (2005); Nesseris and Perivolaropoulos (2005)) 

XsNei^M, W) = A(Qm, w) . (4) 
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where 



B[\Im,w) = 2 ' 

i ^» 

c = Ei^- (5) 

Here subscript i denotes the ith SNe la data and (7j is the observed uncertainty. The fi'-' and 
/i"^ are the observed and theoretical distance moduh, respectively. Note that the /if (^i) is 
independent on Hq and given by 

fj,J{zi)^5\og{HodL{zi)) (6) 



2.2. Bciryon Acoustic Oscillations 

The observations of baryon acoustic oscillations allow to constrain the distance-redshift 
relation at different epochs, which is continuously improved by the Sloan Digital Sky Survey 
(SDSS; York et al. 2000; Eisenstein et al. 2005; Davis et al. 2007). From the final set of 
galaxies observed using the original SDSS target selection criteria together with the 2-degree 
field Galaxy Redshift Survey data (CoUess et al 2003; Cole et al 2005), Percival et al. (2010) 
released the new distance measurements at the redshift z = 0.2 and 0.35. This combined 
sample is composed of 893,319 galaxies over 9100 square degrees in the redshift range 
to 0.5. The BAO relevant distance measure is modeled by the so-called volume distance, 
Dv{z) = [D\{z) z / H {z)Y^^ , with comoving angular diameter distance Da{z) = dz' /H{z'). 
Define the distance ratio = rs{zd)/ Dv{z), where z^ is the redshift of the baryon drag epoch, 
which can be calculated from the fitting formula proposed by Eisenstein and Hu (1998) 

where 

bi = 0.313(n^/i2)-°-^^^[l + 0.607(f]^/i2)°-^^% 62 = 0.238{nmhY^^'' (8) 
The observed distance ratios are given as: 

do.2 = 0.1905 ± 0.0061, 4.35 = 0.1097 ± 0.0036. (9) 
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and the inverse covariance matrix is given by 



Then the statistic is 



^BAO 



30124 
-17227 



-17227 
86977 



(10) 



2.3. Cosmic Microwave Background 

The CMB radiation is the cooled remnant of the hot photon-electron plasma in the 
decoupling epoch of our universe and its precise measurements are critical to the cosmology. 
The WMAP satellite has been measuring CMB temperature and polarization anisotropics 
over the full sky since 2001. After 7 years of observations, Komatsu et al. (2009) provided 
the latest fit results of WMAP Distance Priors, including the acoustic scale I a-, the shift 
parameter i?, and the redshift of the decoupling epoch of photons z^c- The first two param- 
eters capture most of the constraining power of the WMAP data for dark energy properties, 
which are given by 

rs{zdc) 



R{zdc) = ^J^mHlDA^Zdc). (12) 

Here the comoving sound horizon size at the decoupling epoch rs{zdc) is given by (Eisenstein 
and Hu 1998) 

with Qfc and the present-day baryon and photon density parameters, respectively. In this 
paper, we fix = 2.469 x 10~^/i~^ (for T^mh — '2.725K), which is the best-fit values given 
by the 7-year WMAP observations (Komatsu et al. 2009). For the decoupling redshift Zdc, 
we use the fitting formula of Hu and Sugiyama (1996) 

Zdc = 1048[1 + 0.00124(f]b/i2)-o-738][i + g.^Q^h^^^, (14) 

where 

_ oms^iObh^y^ 0.560 

~ 1 + 39.5(^6/12)0.763' 92 - l + 21.l(Qb/l2)1.81- ^ ) 
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Using above WMAP distance priors, the CMB ^ statistic is given by 

xImb = - - ^J) (16) 

where — (Z^, -R, ^dc) are the values predicted by a model and are the corresponding 
maximum-likehhood observe values. The three distance priors obtained from the WMAP 
7-year are given as 

= 302.09 ± 0.76, i?= 1.725 ±0.018, z^c = 1091.3 ± 0.91. (17) 
These results are correlated and the inverse covariance matrix is given by 

/ 2.305 29.698 -1.333 \ 
(C-i) = 29.698 6825.27 -113.180 . (18) 
\ -1.333 -113.180 3.414 / 



2.4. Hubble Constant Eq 

The Hubble constant ifo means the present expansion rate of our universe and is very 
important for determining the cosmological distances. Using the three observations in the 
local universe including the "maser galaxy" NGC 4258, the Cepheid variables and the SNe 
la, Riess et al. (2009) provided a precise measurement result B.^ — 74.2 ± 3.6kms~^ Mpc~^. 
In their analysis, the NGC 4258, located at about 7.2Mpc away, plays as the anchor galaxy, 
for the overall uncertainty of its geometric distance is very small, only 3% by far(Herrnstein 
et al. 1999; Humphreys et al. 2008; Greenhill et al. 2009). On the other hand, the 
240 Cepheid variables obtained with the Hubble Space Telescope (Macri et al. 2006) are 
distributed across six recent hosts of SNe la and the NGC 4258, allowing to directly calibrate 
the peak luminosities of the SNe la, which is crucial to constrain the Hubble constant 
through the SNe la data. The authors had rejected those objects with 0.75 mag above 
the best fitted period-luminosity relation, a error value beyond the "normal" observation 
uncertainty, and the Cepheid sample is reduced to 209. In the SNe la data analysis, only 
140 nearby SNe la with the redshift 0.023 < 2; < 0.1 (Hicken et al. 2009) are used and then, 
utilizing the derivations of the scale factor a, the related luminosity distance is expanded 
to the polynomial function of the redshift keeping highest to order. Finally with the 
well-controlled systematic errors, the SHOES team provided the determination of the Hubble 
constant Hq to ~ 5% precision. Then the statistic is simply given by 

2 [i/o-74.2]^ 
XHubi^o) (19) 
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2.5. Joint Analysis of SNe la, BAO, CMB and Hq Data 

For the four independent observations, the hkehhood function of a joint analysis is just 
given by 

L — I/SNc X LbAO X -^CMB X -^Hub 

= exp[-(xiNe + XBAO+XCMB + XHub)/2]- (20) 

3. Constraints On Time- Varying Dark Energy Models 

In this section, we investigate our two-parameter models of dark energy EOS from the 
joint analysis of (SNe + BAO + BAO + Hq). The observational data depend on w{z) only 
through a complicated integral relation, which effectively smooth out partial differences of 
different models. Consequently, to distinguish different dark energy models from observa- 
tional data is quite difficult. In fact, though the use of more than two parameters can offer 
flexible function form to test the dark energy EOS, the observational constraints on those 
parameters are quite poor. Therefore, in this paper we focus our efforts on two-parameter 
models by carefully analysizing the low {z 0) and high redshift {z oo) behaviors of 
dark energy EOS w{z), while the certain function form determines its evolution with time. 
Our two-parameter models are constructed with considering three limitations: firstly, the 
dark energy EOS w{z) has finite boundary values at 2; = and z — 00; secondly, the w{z) 
function has no more than one local extremum point in the range of redshift z = (0,oo); 
thirdly, the basic bricks of the function are z, 1/(1 -|- ^)(i.c. a) and z/{l + ^)(i.e. (1 — a)). 
In the data analysis, only two parameters are variable, and the possible relevant extra pa- 
rameter is fixed, which enlarges our model space. There is an alternative way to choose two 
parameters: assuming that w{z) can be approximated at low redshift z < 1 and use the first 
two terms of its Taylor series about z — (Copeland et al. 2006). But when considering 
data from the CMB, the high redshift hmit becomes important and our choice is reasonable. 



3.1. Two-Parameter Dark Energy Models 

Though the certain knowledge of w{z) is still lacked and the dark energy modeling 
space is infinite dimensional, our aim is to find some interesting modeling forms and study 
the corresponding behaviors of dark energy. We try different varieties of function forms with 
paying much attentions to keep the models as simple as possible. According to the boundary 
behaviors of the EOS w{z) at 2; = and 2; = 00 and its local extremum property in the range 
of redshift z — (0, 00), we classify eight parameterations of dark energy EOS w{z) into two 
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types of dark energy models, where three of the models have no extremum point and behave 
analogously to the CPL model at two redshift boundaries z — and z = oo. For other five 
models, they all have one local extremum point with the same asymptotic values at the low 
and high redshift limitations. 

Let us first consider three models which are analogous to the CPL model and classified 
as the type I models: 



Just like CPL model, the three Type I models are all monotonic functions of redshift z and 
have the similar boundary behaviors: w{z — 0) — Wq and w{z — >■ oo) = wq + Wi. Thus 
their variety in the whole domain is \wi\, which implies a somewhat slow evolution of dark 
energy. The dark energy EOS w{z) of the type 1 models have no local extremum point and 
its evolution with time is simple. In fact, the cosmological constant or the vacuum energy is 
the simplest model with the fixed EOS value w = —1, and then the constant model with a 
certain value unchanging with time. 

As a simple extension to the type 1 and CPL models, we are now considering models in 
which there exist one local extremum point for the dark energy EOS. In this case, the EOS 
w{z) is no longer monotonic and its derivative cross zero value in the past time. Such models 
are classified as type II models in our present considerations. For the evolution pattern of the 
dark energy, we pay much attention to the logarithmic function forms. As it is known that 
the dark energy EOS of some scalar field models can be approximated by the logarithmic 
function of redshift z, and on the other hand, the quantum corrections of field theory are 
often given by the logarithmic forms. Therefore, the dark energy EOS could be related to 
the logarithmic function of the scale factor a. By using the logarithmic form of a = 1/(1 + ^) 
and (1 — a) = z/{l + z), we can construct the following four kinds of models 



la : w{z) 
lb : w{z) 
Ic : w{z) 
CPL : w{z) 



Wo + Wiz/{1 +z)[l + 1/(1 + z)] 
Wo + ■u;i2;ln[(l + z)/z] 
Wo + wi{l - ln(l + z)/z) 
Wo + Wiz/{1 + z) 



(21) 
(22) 
(23) 
(24) 



lib 



Ild: 



Ila 



lie 



w{z) 
w{z) 
w{z) 
w{z) 



«;o + wi[l/(l + ^)]"ln[l/(l + z)] 

Wo + Wi[z/{l + z)f\n[z/{l + z)] 
wo + wi\n{l + z)/z" 
wo + wiz''ln[{l + z)/z] 



a 



(25) 
(26) 
(27) 
(28) 



where a is an extra parameter and fixed in the data analysis. Note that Model Ic can be 
obtained through the {wq, wi) redefinition of the model lie, and the model lb is the special 
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case of model Ild with a = 1. For a = 0, the model Ila and model lib are reduced to the 
constant model, and the model lie and model lid are reduced to the simple models with 
ln(a) and ln(l — a), which are divergent at boundaries a = and a = 1 or z = oo and z = 0. 
Therefore, the case of a = will not be discussed here. For certain values of a, all models 
have one local extremum in the range 2 = to 2; — ?• oo, and its position is determined by 
the parameter a. The boundary behaviors of the models at high and low redshift depend on 
the parameter a. According to different choices of a, there are three boundary values: Wq, 
Wq + Wi and 00. In this work, we are only interested in the former two cases. For comparison, 
we also present one model without logarithmic function form 

IIe:w;(^) = Wo + Wiz/ {1 + z)" (29) 

where a is an extra parameter and will be fixed in the data analysis. For a = 0, the model 
He is reduced to w{z) = Wq + WiZ, which is divergent when z ^ 00 and will not be considered 
in this work. The corresponding local extremum property is clear that the model He has 
one local extremum point for a > 1. With different choices of a, there are three boundary 
values: wq, wq + wi and 00. For the cases a = 1 and a = 2, it was studied in literature 
(Jassal, Bagla and Padmanabhan 2005). Here we shall consider a as an arbitrary parameter 
to be fixed. In the data analysis late on, we will show that the local extremum property and 
the boundary behavior together determine the pattern of w{z) for its allowed region. 

It is interesting to find in the joint analysis that there exist two strongly correlative 
parameterations 

A : w{z) = Wo + wiz^/^\n[{l + z)/z] (30) 

B:w{z) = Wo + WiZ-^/Hn{l + z), (31) 

which are the special cases of model lid and model lie with a = 4/9 and 5/6, respectively. 
In these two cases, the resulting wq and Wi have almost linear negative correlations with 

Wi = at Wq = —1. 

The boundary behavior at low redshift 2; — )■ and high redshift 2; — )■ 00 for all models are 
summarized in the table [H For the three type I models, the boundary values of dark energy 
EOS w{z) are similar with the GPL model: w{z — 0) — )■ and w{z — )■ 00) — )■ {wq + Wi). 
As type I models are monotonic function, their range is just the inteval wq and {wq + wi). 
For the four models with the logarithmic form, the boundary behavior relies on the value of 
a. For models Ila and lib, their boundary values at the low and high redshift are both wq 
when a > 0. For models lie an lid, their boundary values are both wq for < a < 1. For 
the case of a = 1, the model lie has the low redshift value {wq + wi) and high redshift value 
Wq. For model He, the two boundary values are both wq when a > 1. For the two strongly 
correlative models, their boundary values are also both wq. 
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The local extremum properties of all models have been summarized in Table |2l Though 
the interested model domain in this work is just for the range with 2; > 0, here the redshift 
relating to the future with — 1 < z < is also included and the local extremum point is 
given. For model Ila or model lib, there is one local extremum point located at z* = e^l'^ — 1 
or z = l/(e^/" - 1) for a > 0. Model Ila has one local extremum value in future for a < 0. 
The extremum properties of model lie is somewhat complex. One local extremum point is 
located at z* /{1 + z*) = a\n{l + z*) for < a < 1, and for the case a > 1 one local extremum 
point can be found in the redshift range — 1 < 2; < 0. Model Ild has one local extremum 
value at 1/(1 + z*) = aln[(l + z*)/z*] for < a < 1. For the model He without logarithmic 
function, the local extremum point is found to be at = l/(a; — 1) for a > 1. For the 
two strongly correlative models, their local extremum positions can be obtained from above 
general cases with specific values of a, the corresponding values of the redshift are found to 
be z* = 0.17 for model A and z* = 0.46 for model B. 

The different choice of the extra parameter a will change the boundary values and local 
extremum position of dark energy EOS. For the type II dark energy models, we are interested 
in the cases with finite boundary values and one local extremum point aX z > 0. From the 
summaries given in two tables, the interesting ranges for the extra parameter a are found to 
be: a > for model Ila and mode lib, 1 < a < 1 for model lie and model lid, and a > 1 
for model He. In these cases, the low redshift and high redshift values for all the type II five 
models are equal: w{z — t- 0) = w{z — t- 00) = Wq. 



3.2. Results Prom Latest Observations 

Based on the above analysis, we present in table [3] the best-fit results of five parame- 
ters {wq, wi, Qm, h, lOOflbh^) and the corresponding minimum Xmm' well as the constraint 
results of (wq, wi) with other three cosmological parameters marginalized from the joint anal- 
ysis of (SNe + BAO + CMB + Hq). It can been seen that all Xmm close to Xmm ~ 469, 
which indicates that the current observational data can not explicitly distinguish different 
two-parameters models from the Xmin analysis. For all models, the constraint results for 
three of the cosmological parameters {Qm, h, lOOfi^/i^) are close to (0.28, 0.70, 2.24), though 
the best-fit results for the two parameters wq and wi are quite different. For models Ila, lib, 
and He, the best-fit values of {wo,Wi) vary regularly according to the change of the extra 
parameter a. For models lie and lid, there is a turnaround point in the range 2/3 < z < 1 
and 1/3 < z < 1/2, respectively. 

Let us first discuss the constraints for the three type I models, and compare them with 
the GPL model. In figured], we show the constraint results on {wo,Wi) and w{z) from the 
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OO 
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OO 


Wo 
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w{z) = Wo + Wiz/ (1 + 
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Wq 
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Wq 


Wo + Wi 
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Wo 


OO 


A 
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Table 1: The boundary behaviors of all models 
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z* = 0.46 



Table 2: The local extremum point properties of all models 
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[Wq^ Wi) 
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(-0.97, -0.19), ( 0.28, 0.70, 2.24) 


469.045 


(-o.97l°:?i,-o.i5i;i) 


la 




(-0.97, -0.15), ( 0.28, 0.70, 2.24) 


469.027 


(-0.96l°i?,-0.12;°i^) 


lb 




(-0.95, -0.14), ( 0.28, 0.70, 2.24) 


468.992 


(-0.93l°:i,-0.19;?;ol) 


Ic 




(-0.93, -0.19), ( 0.28, 0.70, 2.24) 


468.959 


(-0.971°:!?, +0.22l°;«^) 
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1 


(-0.9 h 0.33). ( 0.28. 0.70. 2.23) 


1(38.971 


(-0.9r,+"-2T +o.32+}-t') 

V — (J.zo ' — i.ii/ 


2 


(-0.86, 0.63), ( 0.28, 0.70, 2.23) 


468.732 


(-0.86+°-!, +0.64+|-^t) 


3 


(-0.66, 1.33), ( 0.29, 0.69, 2.24) 


468.199 


(-0.66+°f ,+1.32+^-29) 


4 


(-0.57, 1.61), ( 0.29, 0.69, 2.24) 


467.908 


(-0.57+°-^° +1.58+^-^?) 
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1 


(-0.77, 0.87), ( 0.28, 0.70, 2.23) 


468.644 


(-0.77+°-^^ ,0 87+2.28) 


2 


(-0.97, 0.22), ( 0.28, 0.70, 2.24) 


469.023 


(_0.96+°-|L +0.28+?-^^) 


3 


(-0.99, 0.03), ( 0.28, 0.70, 2.24) 


469.077 


(_0.99+°-|2, +0.02+?i2N 


4 


(-1.00, -0.37), ( 0.28, 0.70, 2.24) 


469.074 


(_0.99+°-??, -0.04+?-^3) 

V — U.il' —1.79/ 


TT 

lie 


1/5 


(-0.96, -0.11), ( 0.28, 0.70, 2.24) 


469.026 


(_0.97+o-?i -0.10+?-^!^) 

\ " 0.19' " " — 1.02/ 


1/3 


(-0.95, -0.15), ( 0.28, 0.70, 2.24) 


469.000 


(_0.94+°-^o -0.16+?-^^) 

\ " — 0.26' " " — 1.05/ 


1/2 


(-0.90, -0.23), ( 0.28, 0.70, 2.24) 


468.949 


(_0.91+o-f , -0.22+?-^?) 

V " — O.JY' " — 1.21/ 


2/3 


(-0.79, -0.41), ( 0.28, 0.70, 2.23) 


468.862 


(_0.80+°-^i -0.38+j-^^) 


11(1 


1/5 


(-1.14, 0.12), ( 0.28, 0.70, 2.24) 


469.956 


(-1.13+°i^,+0.10+°i^) 


1/3 


(-1.10, 0.11), ( 0.28, 0.70, 2.24) 


469.053 


(-0.97+?-™, -0.02+J-^^) 


2/3 


(-0.67, -0.58), ( 0.28, 0.70, 2.23) 


468.719 


(-0.681° j^, -0.5611:^^) 


4/5 


(-0.84, -0.33), ( 0.28, 0.70, 2.23) 


468.857 


(-0.87l°f , -0.281?:^?) 


He 


5/4 


(-0.96, -0.22), ( 0.28, 0.70, 2.24) 


469.006 


(_0.96l°it, -0.211}:°^) 


4/3 


(-0.95, -0.25), ( 0.28, 0.70, 2.24) 


468.997 


(-0.951°:^^, -0.241}:^^) 


3/2 


(-0.95, -0.32), ( 0.28, 0.70, 2.24) 


468.977 


(-0.94l°:i,-0.3ll}i^) 


7/4 


(-0.93, -0.45), ( 0.28, 0.70, 2.24) 


468.939 


(-0.93l°:i, -0.441^;^^) 


A 


4/9 


(-0.57, -0.51), ( 0.28, 0.69, 2.24) 


468.947 


(-0.57l°:i, -o.52i?:««) 


B 


5/6 


(0.16, -1.74), ( 0.29, 0.69, 2.24) 


468.525 


(+0.15l°i«,-1.72l?ig) 



Table 3: Based on the joint analysis of (SNe -|- BAO -|- CMB -|- ifo), the best-fit results of 
five parameters (wq, wi, Om, ^, lOOO;,^^) and the corresponding minimum x^j„, as well as 
the constraint results of {wq,Wi) with other three cosmological parameters marginalized. 
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joint analysis of (SNe + BAO + CMB + Hq). The crosshairs in top panels mark the best-fit 
points {wo,Wi) = (-0.96,-0.12), (-0.93,-0.19), (-0.97,-0.22) and (-0.97,-0.15), from 
left to right. In bottom panels, the best fit results are shown as solid lines. Clearly, these 
panels display the following common features. Firstly, it shows negative dependence between 
the parameters wq and wi, which indicates why the area of allowed region for w{z) below 
the best-fit line is larger than the one above the best-fit line. Secondly, around z ~ 0.2, 
the allowed region for w{z) has a narrow node, at which the uncertainty of w{z) becomes 
minimum due to correlations of the parameters Wq and Wi (see, e.g., hinder and Huterer 
2005; Albrecht et al. 2006). Thirdly, the constraining power at high redshift 2; >> 1 is 
weaker than that at low redshift z ~ 1. The reason is that there are more observational 
data at low redshift than at high redshift and the relative low density of dark energy at 
high redshift makes the data insensitive to the model variety. Finally, the resulting best-fit 
w{z) for all the modles is similar with a value around w ~ — 1. It can be seen that the 
constraining power on w{z) is different for various models and the model la considered in 
our present paper leads to the best constraint results. 

We now consider the type II models based on the joint analysis of four observations (SNe 
+ BAO + CMB + Hq). Figure 12] plots the 68% and 95% confidence contours of [wq.Wi) 
and w{z) for the model Ila with the given extra parameter a = 1,2,3 and 4, respectively. 
The crosshairs in top panels mark the best-fit points (wq, wi) = (—0.95, 0.32), (—0.86, 0.64), 
(—0.66,1.32) and (—0.57,1.58), from left to right. The solid lines in bottom panels show 
the best-fit results. It shows the following common features Firstly, the two parameters Wq 
and Wi have positive correlations. Secondly, the confidence intervals of Wi is much bigger 
than that of wo, which means that the observational data are less sensitive to wi. Thirdly, 
there exist two narrow nodes in the allowed regions and the node corresponding to the small 
redshift is near z ~ 0.1. Such a feature is caused by the local extremum point in the model 
Ila around z = 1. Fourthly, the constraining power at low redshift z << 1 is similar to 
that at high redshift z » 1, which is partially due to the same boundary behaviors of low 
redshift and high redshift: w{z — > 0) = w{z — ?■ 00) — ?■ wq- When the parameter a is taken 
to be a larger value, the allowed region between two narrow nodes goes smaller, the best-fit 
line curves down, and the relative constraining power becomes weaker at both low redshift 
and high redshift. 

We plot in figure |3]the 68% and 95% confidence contours for the model lib with various 
values of parameter a = 1,2,3 and 4. The crosshairs in top panels mark the best-fit points 
{wo,Wi) = (-0.77,0.87), (-0.96,0.28), (-0.99,0.02) and (-0.99,-0.04), from left to right. 
The solid lines in bottom panels show the best-fit results. It is not difficult to observe from 
the figure 3 the following common characters: Firstly, the allowed range of wi is much bigger 
than Wo and they have strong positive correlations. Secondly, there are two narrow nodes 
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in the allowed regions. Thirdly, the constraining power at low redshift 2; << 1 is analogous 
to that at high redshift z » 1, which can be explained by the same boundary behaviors 
of low redshift and high redshift: w{z — )■ 0) = w{z 00) — )■ wq. For a larger value of 
parameter a, the allowed interval of w{z) between two narrow nodes increases rapidly, and 
the constraint power gets stronger for low and high redsfhits, which just opposites to the 
case in the model Ila. Especially, the behavior of phantom-like dark energy model indicated 
by the case with a = 4 is consistent with the current observations: the EOS w{z) is close to 
— 1 in both low redshift z < 0.5 and high redshift z > 100 and decreases to much negative 
values when z ^ 4. 

In figure HI we present the allowed regions for the model lie with the values of extra 
parameter a = 1/5, 1/3, 1/2 and 2/3. The crosshairs in top panels mark the best-fit points 
{wo,Wi) = (-0.97,-0.10), (-0.94,-0.16), (-0.91,-0.22) and (-0.80,-0.38), from left to 
right, respectively. The solid lines in bottom panels show the best-fit results of w{z). It 
is different from the models Ila and lib but similar to the type I models that the two 
parameters wo and wi have negative correlations. It can been seen that for a larger value of 
a, the pattern of allowed region changes remarkably, the allowed region between two narrow 
nodes becomes smaller and the allowed regions on both sides of redshift are enlarged. 

The constraint results of the model Ild are shown in figure |5] with values a = 1/5, 1/3, 
2/3 and 4/5. The crosshairs in top panels mark the best-fit points {wq, wi) = (—1.13, —0.10), 
(-0.97,-0.02), (-0.68,-0.56) and (-0.87,-0.28), from left to right. The solid lines in 
bottom panels show the best-fit results of w{z). It is seen that the two parameters Wq and 
Wi have negative correlations. For different values of a, the pattern of allowed region changes 
rapidly and the best-fit values of {wo,wi) have an extremum point in the range < a < 1. 

Figure E] shows the constraint results of the model He with a = 5/4, 4/3, 3/2, 7/4, from 
left to right. The crosshairs in top panels mark the best-fit points {wq, Wi) = (—0.96, —0.21), 
(-0.95, -0.24)), (-0.94, -0.31) and (-0.93, -0.44). The solid lines in bottom panels show 
the best-fit results. It is seen that in all cases there are two minimum values of allowed 
regions and the one with small redshift is fixed to be around z ~ 0.2, and the constraining 
power at low redshift 2; << 1 is similar to that at high redshift z » 1. The allowed patterns 
are relatively insensitive to the change of a. 

3.3. Parameter Strongly Correlative Models 

In two parameter models, it was tried to find the uncorrelated parameters through the 
redefinitions of the two initial parameters (Albrecht et al. 2006; Wang 2008). Alternatively, 
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from the constraint results of above two-parameter models with varying the extra parameter 
a, we find two special cases with very strong correlation between the two parameters Wq and 
Wi- It is shown in figure [7] that two models with strong parameter correlation are found from 
the constraint results to be 

w{z) = Wo + wiz^^^\n[{l + z) / z] (32) 
w{z) = wo + wi\n{l + z)/z^^^. (33) 

In the left panels, we plot the 68% and 95% confidence contours of {wq, Wi) with the crosshairs 
marking the best-fit points {wo,Wi) = (—0.57,0.52) and (0.15, —1.72), from top to bottom. 
It can be seen that the parameters wq and wi are strongly correlated. The linear least-squares 
fitting gives the following relations 

wi = -1.2(1 + Wo) (34) 

wi = -1.5(1 + Wo). (35) 

It is interesting to notice that at Wq = —1, then Wi = and w{z) = —1. From the two right 
panels, it can be seen that in both models, the allowed regions between two narrow nodes 
around z ~ 1 are quite small and the constraining powers at low and high redshift become 
weaker. The situation is in general expected from the following considerations: (a) in the 
region with very low redshift z << 1, the effect of dark energy on the space-time geometry 
can not be accumulated enough to show the differences among different models; (b) at high 
redshift z » 1, due to the relative low density of dark energy based on the observations for 
the evolution of the early universe, it makes the data fairly insensitive to the model variety; 
(c) in the region around 2; ~ 1, there are more observational data and also the different 
dark energy models lead obviously to different cosmological distances. Therefore, the most 
sensitive redshift region of cosmological observations on dark energy EOS w{z) is around 
z ~ 1, which is corresponding to the most strong constraining power. 



4. CONCLUSIONS 

We have investigated a set of time- varying dark energy models by parametering the EOS 
w{z) via two parameters and some of them depend on a given extra parameter. According 
to the boundary behavior and the local extremum point of the EOS w{z), they are classified 
into two types of models. From the joint analysis of four observations (SNe + BAO + CMB 
+ Hq), we have presented the constraint results for all the time- varying dark energy EOS 
w{z). 

Our two-parameter models are constructed with the simple conditions that the dark 
energy EOS w{z) has finite boundary values at z = and z = 00, and the w{z) function has 
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no more than one local extremum point in the redshift range < z < oo. The low redshift 
{z — >■ 0) and high redshift {z oo) behaviors of dark energy EOS w{z) are described by two 
parameters, and the local extremum property and boundary behavior together characterize 
the pattern of w{z) for its allowed region. For different given values of an extra parameter 
a, we have carefully studied the w{z) evolution and found out its reliable values. 

Based on the joint analysis of (SNe + BAO + CMB + Hq), we have shown for all 
models the best-fit results of five parameters {wq, Wi, ^Im, h, lOOQiih"^) and the corresponding 
minimum Xmin^ ^^^1 constraint results of {wo,Wi) with other three cosmological 

parameters marginalized. It have been seen that all Xmin close to Xmiu — 469, which indi- 
cates that the present observational data are not sufficient to distinguish the two-parameters 
models from the Xmin analysis. For all models, the constraint results for three cosmological 
parameters (JIm, h, lOOQfe/i^) are close to (0.28, 0.70, 2.24), though the best-fit results of Wq 
andwi are quite different. 

The three type I models have been shown to have similar functional behaviors with the 
CPL model and the constraint results show their common features: (a) the two parameters 
Wo and Wi have negative correlation; (b) there is a narrow node at z ~ 0.2 in the allowed 
regions; (c) the constraining power at high redshift z — > oo is somewhat weaker than that 
at low redshift z ^ \\ (d) the best-fit wi^z^ is quite fiat and close to —1. It has been seen 
that our model la gets a stronger constraint than the CPL model, thus the Model la may 
provide an interesting dark energy model from the best constraint results. 

For the type II dark energy models which have a local extremum point at ^ > 0, we have 
shown the following interesting features for four of the models with the logarithmic function 
form: (a) the confidence intervals of w\ is much bigger than that of w^^ which indicates that 
the observational data are less sensitive to w^] (b) the allowed regions have two minimum 
points in the most of cases; (c) for the model lla, the fitted wi^z) patterns are not sensitive to 
the change of the extra parameter a; (d) when changing the extra parameter a, the variety 
of the allowed region between two narrow nodes is opposite to the regions on both sides. For 
the model He without involving logarithmic function form, it has been found that the fitted 
w{z) patterns are relatively insensitive to the extra parameter a. 

From the joint analysis, we have found two interesting models which correspond to two 
special values of the extra parameter a in models lie and Ild, where the two parameters 
and W\ get strong correlation with almost a linear relation: w\ oc (l-t-Wo). It is also seen that 
the allowed regions of w{z) between two narrow nodes are quite small and the constraining 
powers at low and high redshift become weak. The situation is understood from the fact that: 
(a) in the very low redshift z «\^ the effect of dark energy on the space-time geometry 
cannot be accumulated enough to show the differences among different models; (b) at high 
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redshift z » 1, the quite low density of dark energy, which is required by the observed 
evolution of the early universe, makes the data fairly insensitive to the model variety; (c) 
the most observational data are available in the region around z ~ 1, where different dark 
energy models can obviously lead to different cosmological distances. 

In conclusion, we have introduced a set of two-parameter dark energy models for the 

EOS w{z) and made a joint analysis based on the four observations (SNe + BAO + BAO + 
Hq). We have shown that the constraint results for all models get almost the same Xmin 
the cosmological parameters (I^Mj ^? although the allowed regions for the EOS w{z) 

are sensitive to different models and a given extra model parameter. In particular, we have 
found two models in which two parameters wq and wi get strong correlation and are given 
almost by a linear relation wi oc (1 + Wq). 
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Fig. 1. — The 68% and 95% confidence contours of four models: w{z) = wq + Wiz/{1 + 
z)[2 - z/{l + z)], w{z) = Wo + wiz\n[{l + z)/z], w{z) = Wq + Wi{l - ln(l + z)/z), w{z) = 
Wq + Wiz/{1 + z), from left to right. The solid lines in bottom panels show the best fit 
results. The crosshairs in top panels mark the best-fit points {wo,wi) = (—0.96,-0.12), 

Q.-^ -n ^Q^ (-O Q7 -n 9.9.) and (-O Q7 -n l.^l fmm Isft. tn rio-ht 
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Fig. 2. — The 68% and 95% confidence contours of the logarithmic form model w{z) = 
Wq + Wi[l/(1 + 2;)]"ln(l/(l + 2;))" with a = 1,2,3 and 4, from left to right. The solid lines 
in bottoms panels show the best fit results. The crosshairs in top panels mark the best-fit 
points (wo,wi) = (-0.95,0.32), (-0.86,0.64), (-0.66,1.32) and (-0.57,1.58), from left to 

ri crht 
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Fig. 3. — The 68% and 95% confidence contours of the logarithmic form model w{z) = 
Wq + Wi[z/{1 + 2;)]"ln(z/(l + 2;))° with a = 1,2,3 and 4, from left to right. The solid lines 
in bottoms panels show the best fit results. The crosshairs in top panels mark the best-fit 
points {wo,wi) = (-0.77,0.87), (-0.96,0.28), (-0.99,0.02) and (-0.99,-0.04), from left to 

ri crht 



- 24 - 



3 



2- 




-31 I I I I I I I 

-2 -1 0-2 -1 -2 -1 0-2 -1 




-3 -2 -1 1 2 3 -3 -2 -1 1 2 3 -3 -2 -1 1 2 3 -3 -2 -1 1 2 3 



Log(z) Log(z) Log(z) Log(z) 



Fig. 4. — The 68% and 95% confidence contours of tlie logaritlimic form model w{z) = 
Wq + Wiln(l + z)/z°' with a = 1/5, 1/3,1/2 and 2/3, from left to right. The solid lines 
in bottoms panels show the best fit results. The crosshairs in top panels mark the best-fit 
points (wo,wi) = (-0.97,-0.10), (-0.94,-0.16), (-0.91,-0.22) and (-0.80,-0.38), from 

Ifift. tn ricrht 
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Fig. 5. — The 68% and 95% confidence contours of tlie logaritlimic form model w{z) = 
Wq + Wi2"ln[(l + z)/z] with a = 1/5, 1/3,2/3 and 4/5, from left to right. The solid lines 
in bottoms panels show the best fit results. The crosshairs in top panels mark the best-fit 
points (wo,wi) = (-1.13,-0.10), (-0.97,-0.02), (-0.68,-0.56) and (-0.87,-0.28), from 

Ifift. tn ricrht 
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Fig. 6. — The 68% and 95% confidence contours of tlie model w{z) = wq + wiz/{1 + z)"" witli 
a = 5/4, 4/3, 3/2, 7/4, from left to right. The solid lines in bottoms panels show the best 
fit results. The crosshairs in top panels mark the best-fit points (wq, Wi) = (—0.96, —0.21), 
(-0.95, -0.24)), (-0.94, -0.31) and (-0.93, -0.44), from left to right. 



-27- 



w(z) = + z"Log[(1+z)/z]; a=4/9 




Fig. 7. — The 68% and 95% confidence contours of tlie two most correlation models, ^(2;) = 
Wq + Wi2;^/^ln[(l + z)/z] and w{z) = Wq + Wiln(l + z)/z^^^, from top to bottom. The sohd 
lines in the two left panels are the linear least-squares regression line: Wi = — 1.2(1 + wq) in 
top-left panel and wi = —1.5(1 + wo) in bottom-left panel. 
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